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We study the coherence trapping of a qubit correlated initially with a non-Markovian bath in 
a pure dephasing channel. By considering the initial qubit-bath correlation and the bath spectral 
density, we find that the initial qubit-bath correlation can lead to a more efficient coherence trapping 
than that of the initially separable qubit-bath state. The stationary coherence in the long time limit 
can be maximized by optimizing the parameters of the initially correlated qubit-bath state and the 
bath spectral density. In addition, the effects of this initial correlation on the maximal evolution 
speed for the qubit trapped to its stationary coherence state are also explored. 
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Introduction. Quantum state takes the form of su¬ 
perposition which leads to quantum coherence. Quan¬ 
tum coherence plays a central role in the applications of 
quantum physics and quantum information science Bi- 
However, it is fragile due to interactions of the environ¬ 
ment. Understanding of quantum coherence dynamics of 
an open system is a very important task in many areas of 
physics ranging from quantum optics to quantum infor¬ 
mation processing. It is known that many quantum open 
systems exhibit non-Markovian behavior with a flow of 
information from the environment back to the system 
This presence of non-Markovian effects can in¬ 
duce the long-lasting coherence in biological surroundings 
and the steady state entanglement in coherently coupled 
dimer systems ii- By considering the pure dephas¬ 
ing non-Markovian bath, decay of quantum coherence of 
the system would be terminated in a finite time, such 
that the system can partly retain coherence in the long 
time limit. This new phenomenon, known as coherence 
trapping [Toj| , is important for quantum information pro¬ 
cessing since the effective long-time quantum coherence 
of the system is preserved. Coherence trapping of a quan¬ 
tum system is mainly related to the open dynamics, and 
is generally analyzed in the fact that the system and bath 
are initially separable. As is well known, however, the 
initial system-bath correlations are important for the the 
dynamics of the open systems. The distinguishability of 
quantum states would increase in the presence of initial 
system-bath correlations [13, Si • The information flow 
between the system and its bath and the correspond¬ 
ing degree of non-Markovia nity can also be influenced 
by these initial correlations |19l - [^ . On the other hand, 
the standard master equation approach to open systems 
may not be appropriate unless a product state is explic- 
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itly prepared [ 11113 . So the coherence trapping of an 
open system due to the initial system-bath correlations 
should be studied both physically and methodologically. 

In this paper, we will concentrate on the following 
questions: how do the initial system-bath correlations 
affect coherence trapping of the system? which form of 
the initially correlated system-bath state can maximize 
the stationary coherence of the system? We consider the 
pure dephasing model of a qubit initially correlated with 
a zero-temperature Ohmic-Like bath. We will show that 
the initial qubit-bath correlation can lead to the more 
efficient coherence trapping, while the lower initial co¬ 
herence of the qubit is induced by this initial correlation. 
In the long time limit, the stationary coherence of the 
qubit can be maximized by choosing the optimal param¬ 
eters of the initially correlated qubit-bath state and the 
optimal Ohmicity parameter of the bath. 

Furthermore, the task to drive an initial state to a pre¬ 
scribed target state in the shortest possible time is signif¬ 
icant for quantum control in many areas of physics, such 
as quantum com put ation [3 , f^ist population transfer in 
quantum optics [24| . and quantum optimal control pro¬ 
tocols [ 13 , 113 ] ■ This minimum evolution time, which is 
defined as quantum speed limit (QSL) time [27 h 3 - is a 
key method in characterizing the maximal speed of evo¬ 
lution of quantum systems. Here in order to speed up the 
evolution from an initial coherence state to its stationary 
coherence state, we further focus on the interactions of 
the initial qubit-bath correlated state, the spectral den¬ 
sity function of the bath and the QSL time. Remarkably, 
we find that the initial qubit-bath correlation can reduce 
the QSL time for the occurrence of coherence trapping. 
The maximal evolution speed for the qubit trapped to its 
stationary coherence state can also be controlled by opti¬ 
mizing the parameters of the initial qubit-bath correlated 
state and the bath spectral density function. 

Model. Let us consider an exactly solvable model, in 
which the process of energy dissipation is negligible and 
only pure depahsing is a mechanism for decoherence of 
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the qubit. The associated Hamiltonian reads (setting 


H = woCTz + / Loalaujduj + / cr^[g^al + glaj\duill) 
Jo Jo 

where the operator cr^ is defined by (Tz = |e)(e| — \g){g\, 
associated with the upper level |e) and the lower level \g) 
of the qubit; aui and are the bosonic annihilation and 
creation operators for the bath, which is characterized 
by the frequency w; g^i is the coupling constant of the 
interactions of the qubit with the bath, and g* is the 
complex conjugate to The Hamiltonian in Eq. (IT]) 
can be rewritten in the block-diagonal structure [^, 

H = diag[He, Hg], where Hf./g = iwo + /g°° ± 

Here, we consider the situation where a correlated ini¬ 
tial state of the qubit-bath system in the form (l8j |. 


r(t)], where. 
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where r[-] is the Euler gamma function and the parame¬ 
ter i} = lfi + v)/2. 

Coherence trapping for the quhit. How to quantify 
quantum coherence of a quantum system now becomes 
paramountly important. In recent years, a wide vari¬ 
ety of measures of coherence have been proposed . 

Currently, Baumgratz, Cramer and Plenio find that the 
relative entropy of coherence [d^ . 


C{p) = Sipd^ag) - S{p), (4) 


|4'(0)) = Ce|e) O l^o) + Cg\g) O |6), (2) 

with the non-zero complex numbers c^/e Eire satisfied 
|cep + |cgp = 1. And we assume that |^o) is a bath 
ground state and |^a) = — A)|^o) + M^f)] is a 

bath superposition state of the ground state |^o) and a 
coherent state |^/) = D{f)\fo)- The displacement oper¬ 
ator D{f) reads D{f) = exp{/“[/,^ai - /><^]dw} for 
an arbitrary square-integrable function /. The constant 
Cx = x/{l — A)^ -I- -I- 2A(1 — A)i?e(^o|^/) normalizes 

the state |^a), where Re is a real part of (Co 10) i^ i'ii® 
bath Hilbert space. The correlation parameter A G [0,1] 
determines the initial correlation of the qubit and bath. 
Through performing the Hamiltonian described in Eq. 
o, the state of the total system at any time t is given 
by |4'(t)) = Ce|e) O 10(0) + Cglff) C) |0(^))> where 
|0(i)) = exp{-iHet)\^o) and |O(0) = expl-iHgt)\^x)- 
Then the reduced density matrix px{t) of the qubit at 
time t reads, Pee{t) = |cep, Pggit) = |cgp and Peg{t) = 
p*ge{t) = CeC*gTx{t), with the dephasing rate Tx{t). 

The qubit dynamics is closely dictated by the spec¬ 
tral density function characterising the qubit-bath inter¬ 
action. In the following the bath can be described by the 
family of Ohmic-Like spectra \gu}\^ = exp(—w/wc), 

with LOc being the cutoff frequency and a > 0 a dimen¬ 
sionless coupling constant. By changing the /i-parameter, 
one goes from sub-Ohmic baths (—1 < ^ < 0) to Ohmic 
(/r = 0) and super-Ohmic (p > 0) baths, respectively. 
Furthermore, the coherent state |C/) can be calculated by 
the spectral density function = 0 ;®^+^ exp(—w/wc), 
with V > 0. So the initial state of the qubit-bath 
system can be controlled by the parameters A and v. 
For the case A = 0 the qubit and the bath are ini¬ 
tially uncorrelated, the dephasing rate can be obtained, 
To(t) = exp[—2ia;ot—r(t)]. While for 0 < A < 1 the initial 
correlation exists in the qubit-bath system, one also finds, 
Tx{t) = —A-|-Aexp[—2i(;ii(t)-|-fc(t)]}exp[—2*a;ot — 


where S{p) is the von Neumann entropy and pdiag de¬ 
notes the state obtained from p by deleting all off- 
diagonal elements, and the intuitive h norm of coherence, 

C/i(p)= (5) 

i,3J¥=3 

are both general and proper measures of coherence. So 
in the following we will choose the relative entropy of 
coherence C (p) to measure the quantum coherence of the 
reduced density matrix p\{t) of the qubit in the presence 
of qubit-bath initial correlation. 

If there is no correlations in the initial qubit-bath state, 
the qubit dephasing To(t) is characterized by exponential 
decay of the qubit coherence, hence will predict vanishing 
coherence in the long time limit in the Ohmic and sub- 
Ohmic dephasing baths [l^. While for the super-Ohmic 
baths, the qubit dephasing will stop after a finite time, 
therefore leading to coherence trapping. This behavior 
can realize the effective long-time coherence protection. 
In the following, we would mainly see the effect of the 
initially correlated qubit-bath state on coherence trap¬ 
ping of a qubit in the super-Ohmic bath model. The 
preparation of this initially correlated qubit-bath state 
can be obtained by non-local operations with two steps 
p^ . Firstly, one prepares the bath state |^a), and then 
products it with the ground state of the qubit Ip). Then, 
one would superpose |p) O I'Ca) with the product state 
|e) ® 1 ^ 0 )) with the weights Cg/ei respectively. And the 
initial correlations of the qubit-bath system can be con¬ 
trolled by the parameters Cg/e, A and the function /^. 

We shall examine the decoherence process where the 
initially correlated qubit-bath state is in the form of 
Eq. ([2]), with Ce = Cg = Then the initial co¬ 

herence of the qubit can be evaluated C(p^'"°) = ^(1 — 
TA(0))log2[l-TA(0)] + i(l + TA(0))log2[l + TA(0)],with 
TrA(O) = C'a ^(1 — A -I- Aexp[—ir[r;]w)(]). At time t = 0, 
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FIG. 1: (Color online) The stationary coherence of the qubit 
quantified by the relative entropy of coherence C'(poo) as a 
function of the bath parameters a and fi. (a) for the uncor¬ 
related initial qubit-bath state (A = 0); (b) for the correlated 
initial qubit-bath state (A = 1); (c) a = 0.2; (d) A = 1. 
Parameters are chosen as, v = 1.5, and ujc = 1. 


in the case A = 0 the dephasing rate To(0) = 1, while for 
the correlated initial state we can obtain 0 < Ta( 0) < 1. 
So it means ) < C{pf^Q), that is to say the initial 

correlation of the qubit-bath system can lead to lower 
initial coherence of the qubit. 

On the other hand, to clear the effect of the qubit-bath 
initial correlation explicitly, we also perform the calcula¬ 
tion for the stationary value of coherence trapping in the 
long time limit. In Fig. 1, we show the stationary co¬ 
herence C{poo) between the initially uncorrelated A = 0 
and correlated A = 1 states as a function of the bath 
parameters a and p. By comparing Figs. 1(a) and 1(&), 
it is clear that the presence of the qubit-bath correla¬ 
tion in the initial state enlarges the region for occurrence 
of coherence trapping. Moreover, by giving the other 
parameters, Fig. 1(c) clearly shows that the larger cor¬ 
relation parameter A leads to a more efficient coherence 
trapping as the stationary coherence is higher than that 
obtained from the initially uncorrelated qubit-bath state. 
Although the lower initial coherence of the qubit can be 
induced by the correlation parameter A, the coherences of 
the bath subsystem and the qubit-bath composite system 
would appear in the initial qubit-bath state correspond¬ 
ingly. That is the main physical reason of the more effi¬ 
cient coherence trapping of the qubit induced by the cor¬ 
related initial qubit-bath state. Additionally, from Fig. 
1 we also can easily find that, the stronger coupling a of 
the qubit to bath diminishes the stationary coherence in 
the long time limit. And there exists an optimal value of 
the Ohmicity parameter p = 1.46 of the bath maximiz¬ 


ing the stationary coherence, which is independent of the 
coupling constant a and the correlation parameter A, as 
shown in Figs. 1(c) and l{d). 

Next, by choosing the optimal value p = 1.46 of the 
super-Ohmic bath, the influence of the parameters char¬ 
acterizing the initially correlated state on coherence trap¬ 
ping is depicted in Fig. 2(a). Two regions, the enhancing 
of coherence tapping (ECT) (i.e. C(p^°) > C{p^^) = 
0.1827 ) and the no-enhancing of coherence trapping 
(No-ECT) (i.e. C{p^^)<C{p^°)), are acquired in the 
corresponding parameter planes. The dashed-white line 
C(Poo) = 0.1827 is the dividing line between these two 
regions. That is to say, not all but specific initial states 
I^a) can lead to the enhancing coherence trapping. The 
range of v to gain the enhancing of coherence trapping, 
would reduce as the correlation parameter A increasing, 
as shown in Fig. 2{b). So we conclude that, in order 
to achieve the most efficient coherence in the long time 
limit, both the optimal Ohmicity parameter p and the 
optimal state |^a) must be satisfied. 



FIG. 2: (Golor online) The stationary coherence of the qubit 
quantified by the relative entropy of coherence C{poo) as a 
function of the parameters for the initial qubit-bath state A 
and V. The dashed-white line in (a) means C(poo) = 0.1827, 
which is the dividing line between two regions. Parameters 
are chosen as, a = 0.2, p = 1.46, and ujc = 1. 


Quantum evolution speed. Since the qubit would occur 
coherence trapping after a hnite time tc in the super- 
Ohmic bath, one may naturally concern the evolution 
speed between the initial state p\{0) and the stationary 
coherence state p\{tc). The quantum speed of evolution 
from pa(0 ) to its target state p\{tc) can be characterized 
by QSL time The definition of QSL time between 

an arbitrary initially mixed state po and its target state 
Pr, governed by the master equation pt = Ltpt, with Lt 
the positive generator of the dynamical semigroup, is as 
follows [1^ tqsl = max{—i , ^ j B{po,pr), 

with A = T-^fJXdt, B(po,Pr) = Itr(poPr) - tr{pl)\ 
denotes a metric on the space of the initial state po and 
the target state pr via the so-called relative purity, and 
ai are the singular values of pt and pi those of the initial 
mixed state po- The above expression of tqsl can ef¬ 
fectually define the minimal evolution time for arbitrary 
initial states, and also be used to assess quantum evolu¬ 
tion speed of open quantum system. 
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FIG. 3: (Color online) The QSL time for the evolution from 
the initial state Pa(0) to the stationary coherence state p\{tc), 
quantified by rqsL/tc as a function of the parameters p and 
V, with ujc = 1. Parameters are chosen as, (a) a = 0.2, v = 2; 
(b) a = 0.2, p = 1.46; (c) a = 0.2, A = 0.5. 

Here, we also consider the weights Ce = Cg = l/-\/2 
in the initially correlated qubit-bath state in Eq. (El- 
Then the QSL time for the qubit initial state px{0) to 
the stationary coherence state p\{tc), can be calculated 

tqsl/U = |Ta(0)[Ta(4)-Ta( 0)]|//J‘= |tA(t)Mt. In Figs. 
3(o) and 3(6), we demonstrate how the QSL time for 
evolution from px(0) to pxitc) can depend on the pa¬ 
rameters p and V, with different selected correlation pa¬ 
rameter A. Firstly, it is clear that the initial qubit-bath 
correlation can reduce the QSL time as the value of A 
increasing. That is to say, the evolution from the initial 
coherence state to the stationary coherence state, can be 
speeded up by the initial correlation in the qubit-bath 
state. And then, another remarkable feature can be ac¬ 
quired: There exist the optimal Ohmicity parameter p 
or the parameter v of |^/), which can induce the mini¬ 
mum value of QSL time. And the optimal parameters 
p OT V are dependent of the correlation parameter A. In 
Fig. 3(a), when v = 2, the optimal Ohmicity parameter 
p = 2.84,2.60,1.80,2.09 for A = 0, 0.3,0.6,0.9, respec¬ 
tively. By choosing p = 1.46, as shown in Fig. 3(6), the 
optimal parameter for the initial bath state |^/) can be 
obtain v = 2.61,2.80,3.09,3.97 for A = 0.1, 0.3,0.6,0.9, 
respectively. 

Furthermore, since both the Ohmicity parameter p and 
the parameter v of |^/) can bring about the minimum 


QSL time, in the following we would seek the optimal 
condition {v,p) on the maximal evolution speed of the 
qubit. Fig. 3(c) shows QSL time for pa( 0) to px{tc) as 
a function of p and v. By a given correlation parame¬ 
ter A = 0.5, we observe that, the minimum QSL time 
can only appear in the region (v < 5,/r < 4). And the 
optimal values {v = 3.65, p, = 3.10) which lead to the 
minimum QSL time Tq^/Ic = 0.2591, can be found by 
accurate numerical calculation. This can be understand 
that, in order to speed up the evolution speed of the 
qubit, the Ohmicity parameter p and the parameter v 
of 1^/) should be optimized. Combined with the above 
section about coherence trapping, the aim to make the 
qubit trap in a higher stationary coherence state with 
the maximal evolution speed, can be attained by choos¬ 
ing the optimal parameters of the initial qubit-bath state 
(A,t>) and the bath spectral density function [p). 

Conclusion. In summary, we studied intriguing fea¬ 
tures of coherence trapping of a qubit with a zero- 
temperature structured bath by considering the initial 
qubit-bath correlation. The initial qubit-bath correla¬ 
tion not only leads to a more efficient coherence trap¬ 
ping, but also speeds up the evolution for the occur¬ 
rence of coherence trapping. Moreover, both the max¬ 
imum stationary coherence in the long time limit and 
the minimum QSL time from the initial state to the 
stationary coherence state, can be acquired by optimiz¬ 
ing the parameters of the initially correlated qubit-bath 
state and the bath spectral density. This physical mech¬ 
anism leading quickly to a higher stationary coherence 
would play an important role for implementing quantum 
simulators [i^ and quantum information processors [i^ . 
It is worth pointing out that the non-Markovian effects 
may not monotonically cause the acceleration of the sys¬ 
tem evolution in the super-Ohmic bath, as shown in Fig. 
3(a). This is clearly different from the main result in the 
damped Jaynes-Cummings model [s^, which shows that 
the evolution speed can be monotonically increased by 
non-Markovian effects. So the specific interplay between 
the evolution speed of the system and the bath non- 
Markovian effects should be studied under different cir¬ 
cumstances. Experimentally, the coherence trapping can 
be demonstrated by qubit-bath systems like optics M , 
trapped ions [i^ and superconducting qubit [13, IH, . 
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